####  Figure S17: Weight distribution ####

# Libraries
# library(tidyverse)
# library(rio)
# library(here)

# data_pnas = import(here("Data","data_pnas.rds"))
# source(here("Code","Functions","funcs_theme.R"))

fig_s17 = data_pnas |> 
  ggplot(aes(x = weight*wt_adj)) + # Use weekly weights
  geom_histogram() +
  theme_prl() +
  labs(x = "Final Weight", y = "Count")

print(fig_s17)

ggsave(here("Plots","Supplementary","figure_s17.png"),
       fig_s17,
       units = "in",
       width = 5, height = 3.5,
       dpi = 600)